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Abstract. A different formulation of the effective interaction hyperspherical harmonics 
(EIHH) method, suitable for non-local potentials, is presented. The EIHH method for local 
interactions is first shortly reviewed to point out the problems of an extension to non-local 
potentials. A viable solution is proposed and, as an application, results on the ground-state 
properties of 4- and 6-nucleon systems are presented. One finds a substantial acceleration in 
the convergence rate of the hyperspherical harmonics series. Perspectives for an application to 
scattering cross sections, via the Lorentz transform method are discussed. 



1. Introduction 

In solving the Schrddinger equation for the ground state of a many body system, via expansion 
techniques, one very often encounters the problem of convergence. In order to deal with this 
problem already several decades ago the notion of effective interaction (EI) has been introduced. 
The true EI is the operator that replaces the bare interaction when working in a finite model 
space, to give (at least) the same ground-state energy. In general it is an ^4-body operator. 
Therefore in the following we will denote it by V^ e ff and refer to it as to the total effective 
interaction. Finding V^ e ff, however, is as difficult as solving the original problem, therefore 
over the years one has searched - more or less successfully - for the best approximations to it. 

A different point of view has been introduced in the nineties with the No Core Shell Model 
(NCSM) PQ. The basic philosophy behind the NCSM success is that the EI is a partial effective 
interaction, tailored for the model space, however, constrained to coincide with the bare one, 
when enlarging the model space. This ensures that the energy will converge faster to the correct 
value. So, instead of pointing at getting the best approximation to the binding energy at a fixed 
model space, one points at reaching the converged solution as fast as possible, increasing the 
model space. 

The same idea is behind the so called EIHH method [2]. The difference is that, instead of 
the harmonic oscillator basis (HO), the hyperspherical harmonics (HH) basis is used. There are 
some advantages in the EIHH approach: the EI is a function of the hyperradius parameter (a 
kind of "hyperlocal" effective interaction) and in addition it is state dependent, resulting in a 
faster convergence, in the relevant quantum numbers, than in NCSM. The disadvantage is that 
working with the HH basis one does not take advantage of the well known versatility of the HO 
basis. In particular the straightforward application of the EIHH method to non-local potential 



presents some problems. In the following it will be shown what these problems are and how 
they can be overcome. 

2. Short summary of the EIHH method for local potentials 

2.1. Hyperspherical coordinates and hyper spherical harmonics 

One characteristic of all the HH methods is that they respect translational invariance. In 
fact the hyperspherical coordinates are defined by a transformation of the Jacobi coordinates 
ffi (i = 1, 2, ... , A — 1), in analogy with the spherical coordinates (r, CI2) of the 2-particle case. 

Of the (3A — 3) coordinates the transformation retains the (2A — 2) Jacobi vector angles. 
The remaining (A — 1) coordinates are obtained generalizing the 3-dimensional space into a 
(A — l)-dimensional space, spanned by the modulus of the (A — 1) Jacobi vectors. Therefore, 
instead of the square radius r 2 , that is the sum of the three squared coordinates, one has a 
square %per-radius, p 2 that is the sum of the square modulus of the (^4 — 1) Jacobi vectors. The 
remaining (^4 — 2) coordinates are again angles and are the generalization of the definition of VL2 
in the 3-dimensional space to &A-2 m the (^4 — l)-dimensional space. Summarizing, the hyper- 
spherical coordinates include one Ziyper-radius p and (3 A — 4) hyper- angles. We will denote all 
hyper-angles by Qa- So any function of the Jacobi coordinates /(jfi, T72, ■ ■■VA-i), when expressed 
in %per-spherical coordinates becomes f(p, &a)- 

The nice feature of these coordinates is that the kinetic energy operator of a many particle 
system takes a form in perfect analogy to the 2-particle 3-dimensional case: a p-dependent 
Laplacian and a /iyper-centrifugal barrier, with a /ij/per-angular momentum operator K that 
depends on all the Ziyper-angles. Therefore the translational invariant Hamiltonian for an A- 
particle system reads. 

^ = -5^ + s;f + ^<*^. m 

The hyperspherical harmonics y\K] are eigenfunctions of K with eigenvalues K{K + 3^4 — 5). 
They constitute a useful basis where one can expand the ^4-particle wave function. It is as well 
useful to expand the hyperradial part in Laguerre polynomials L%(p) so that one has 

*(p,jU) = E c n[K] K( P )y [K] (h A ) , (2) 

[K]n 

Of course the wave function must be complemented by the spin-isospin parts. The whole function 
must be antisymmetric. Therefore one needs to classify the hyperspherical harmonics according 
to the irreducible representations of the symmetry group of A particles. This is a non trivial 
task, that, however, has been solved in [3]. 

2.2. Determination of an effective interaction for local potentials 

In order to construct an EI we use the Lee-Suzuki similarity transformation [4J. This requires 
the definition of an A-particle model space P (Q denotes the rest of the whole Hilbert space, 
P + Q = I). For us P will be spanned by all the A-body HH with K < K max . In order to get 
the total EI the Lee-Suzuki method gives the recipe to get the similarity transformation that 
defines H^ e ff, and therefore V^ e ^ . However, as already said above, we do not search for 
the total EI, but for a partial EI that is tailored for our HH model space and constrained to 
coincide with the bare one, only when enlarging the model space. This is naturally achieved in 
the following way: 

• consider p as a parameter; 



focus on a part of H^ A \ that we will call H®(p), 



H [2] (p) = ^^ + V AA . 1 (p,h A , A ^). (3) 



The term Va,a~i is just any of the pair potential terms in the Hamiltonian. For simplicity 

we choose the pair defining the "simplest' Jacobi coordinate i]A,A-i = \J\{?A — 'f f A-i)-, since 

anyway we are working with antisymmetrized states. The eigenfunctions of H^ 2 \p) will live 
in a space that is contained in the A-body Hilbert space. There will be a P 2 contained in 
P and a corresponding Q 2 contained in Q. Both of them are known since the eigenvalue 
problem for H^ 2 \p) can be easily solved, being a sort of 2-body problem "immersed in a 
medium" (H^ depends on p!); 

• apply the Lee-Suzuki similarity transformation to get the information residing in Q 2 into 
P 2 and obtain the EI V^ e ^ (p); 

• use this EI in the P space; 

• increase the P space to convergence. 

3. The effective interaction for non-local potentials 

For non-local interactions one cannot proceed in the same way as described above. The reason 
is that for non-local interactions Va,a-i is a function of (p, £Ia,A-i', p\ Q'a A-i)- as before the 
variable p is considered a parameter, for fixed p and p' one has Va,a-i(p, &a,A-i', p', &'a A-i) ^ 
Va,a-i(p, Q'a A-v p' 1 &A,A—i)- This means that the potential is not hermitian! 

A solution to overcome this problem is to start again from the full , which for a non-local 
potential is written as 

H,A] = -i A >+i1 ^ M o>.<WA) 

1 IK 2 

= "^ A p + + Vu(P> n ^,p', fiia) + ...V a ,a-i(p, n A ,A-i;p', Ma,a-i) > ( 4 ) 

and subtract and add a term that is local only in p, i.e. Va,a-i(p, &a,a-i', P-, &'a A-i)- 



1 a 1 K 2 

2m p 2m p 2 

1 IK 2 



H^= - ±A p + ^ + vW(p,CiA-,p>,n> A 



= ~ 2m A P + 2 m *p2 +^i2(P>^i2;p / ,^i 2 ) + -Va,a-i{p,&a,a-i;p',&a,a-i) ~ 

- V 12 (p, fii2| P, + V 12(P, ^125 P, ^12) - ••• + ••• 

- V a ,a-i(p,&a,a-i;p,&a,a-i) + V a ,a-i(p,&a,a-i;p,&a,a-i) ■ ( 5 ) 

From this expression one can single out a new (p) formed by the hypercentrifugal term and 
the last term 

H [2] (P) = +V A ,A-l(p,n A ,A-i;p,Cl' A , A -l) ■ (6) 

In this quasi two-body Hamiltonian all non-local effects regarding the hyperspherical coordinates 
are incorporated, while the hyperradial part of the non-locality remains excluded. 



Table 1. Convergence of the HH expansion for the 4 He ground-state energy (in MeV) and root- 
mean-square radius (in fm) with the bare and the effective non-local Idaho N 3 LO potential. 
Results with other methods are also given for comparison. 





Bare 


Effective 


K 

-■'max 


< H > 


V< r 2 > 


< H > 


V< r 2 > 


2 


-3.507 


1.935 


-17.773 


1.620 


4 


-13.356 


1.523 


-22.188 


1.533 


6 


-20.135 


1.446 


-24.228 


1.496 


8 


-23.721 


1.451 


-25.445 


1.498 


10 


-24.617 


1.470 


-25.363 


1.506 


12 


-25.115 


1.491 


-25.439 


1.515 


14 


-25.259 


1.501 


-25.398 


1.516 


16 


-25.310 


1.509 


-25.390 


1.518 


18 


-25.359 


1.513 


-25.385 


1.518 


20 


-25.370 


1.515 


-25.381 


1.518 




-25.37(2) 


1.515(4) 


-25.38(1) 


1.518(1) 


HH [6] 


-25.38 


1.516 






FY [HE] 


-25.37 








NCSM [8] 


-25.39(1) 


1.515(2) 







Applying the Lee-Suzuki transformation to this new HM(p) one obtains VjSf^(p). Then 
the effective non-local interaction is obtained by 

v[ [K*UK> A ] fa P') = V fL],[K' A ) 0>, P) + *(P ~ P) ( V K[*a P) ~ V tLuK> A ] (P. P)) ■ ( 7 ) 

The procedure described here is perfectly in line with the point of view stated above: all what 
is left out in the definition of the partial EI is recovered by increasing the model space P up to 
convergence. In particular, since Eq. (|7|) represents a non-local HH effective interaction derived 
from the diagonal hyper-radial matrix element of in the position representation, it can be 
generalized to an arbitrary hyperradial basis. The question is which of this choices will lead to 
a better effective interaction, i.e. to a faster convergence of the HH expansion. This point has 
been studied in [5]. 

4. A test of the EIHH for non-local interaction 

Here we present applications of the procedure described in the previous chapter to the case 
of 4- and 6-nucleon systems. In Table 1 we present results for ground-state energy and root- 
mean-square radius of 4 He, obtained with the bare non-local Idaho N3LO potential [9] and 
with the corresponding effective interaction, using the hyperradial basis that leads to the best 
convergence [5]. (In order to work with this force we use a representation on a HO basis) 

One sees that Kmax = 8 for ground-state energy and Kmax = 12 for the radius are already 
sufficient to obtain a convergence accuracy of less than 1%. 

For A = 6 nuclei the situation is illustrated in Fig. 1, for the ground-state energies of 6 He 
and 6 Li with the bare JISP16 nuclear force |10] and with the corresponding non-local EI. One 
notices that while it is not possible to obtain converged HH results for the bare interaction, the 
EI convergence is reached with a rather small effective interaction model space. 




5. Perspectives for reaction cross sections 

Calculating a reaction cross section ab initio is considered a much more difficult task than 
calculating binding energies or bound state observables in general. The reason is that in most 
cases reactions involve scattering states. The many-body scattering problem may lack a viable 
solution already for a very small number of constituents in the system. 

The Lorentz Integral Transform method [12] is able to overcome this longstanding 
stumbling block. In fact this approach reduces the scattering problem to a bound state-like 
problem. The essence of the method lays in finding the solution of a Schrodinger-like equation 
with a source S that depends on the reaction one is study. In practice one has to solve 

(H-E -a)\*) = S, (8) 

for many values of the real part <jr of the parameter a = <jr + icrj , and for a fixed value of its 
imaginary part 07, rigorously different from zero. It is just this last condition that ensures that 
the solution of Eq. ([5]) has bound-state boundary conditions. Therefore it can be found by any 
method suitable for bound states, like that described in this contribution. 

The results found in the application of the present version of the EIHH method opens up a 
concrete possibility to study structure and reaction of a few/many-particle system within the 
same framework, also using potentials that are non-local, like the recent effective field theory 
potentials. This allows a critical review of different potential models, whose performances can 
be judged on a much larger set of observables. 
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